COMPUTATION OF SELBERG ZETA FUNCTIONS ON HECKE TRIANGLE 

GROUPS 



FREDRIK STROMBERG 



Abstract. In this paper, a heuristic method to compute the Selberg zeta function for 
Hecke triangle groups, G,j is described. The algorithm is based on the transfer operator 
method and an overview of the relevant background is given. We give numerical support 
for the claim that the method works and can be used to compute the Selberg Zeta func- 
tion on Gq to any desired precision. We also present some numerical results obtained by 
implementing the algorithm. 
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1. Introduction 

The Selberg zeta function, Zp {s), for a co-finite Fuchsian group Y plays an important 
role in the spectral theory or harmonic analysis on the corresponding orbifold ^ = r\J^, 
a surface with constant negative curvature. Selberg's [27J motivation to introduce Zp (i) 
was the similarity between a trace formula he developed (cf. in particular lIZTl p. 74]), now 
called the Selberg trace formula and Weil's explicit formula (31). The role of the Riemann 
zeta function ^ (s) in the latter is analogous to the role of Zp (s) in the former. For a 
more detailed account of this motivation see Hejhal lITOll (in particular sections 4-6). Since 
then the Selberg trace formula has been worked out in detail for PSLai^) (by e.g. Hejhal 
ifTTlfTH ) and the properties of Zp (s) has been extensively studied in many other contexts. 

Despite the importance of Zp (s) and the fact that one can obtain an abundance of its 
properties through the Selberg trace formula, numerical studies of its behavior inside the 
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critical strip, \3is\ < 5 have been surprisingly scarce in the literature. The main reason is 
of course the fact that the defining formula does not represent an analytic function in this 
domain so one is forced to, one way or the other, analytically continue this expression. 

To the authors knowledge, even for the simple case of the modular surface, the only 
successful numerical evaluation of Zp on the critical line was made by Matthies and Steiner 
They overcame the difficulty by desymmetrizing the modular surface Jf /PSLxCE) 
with respect to reflection in the imaginary axis and then restricting their analysis to the odd 
part, which conveniently avoids any interference by the continuous part of the spectra. For 
this system, corresponding to a billiard with Dirichlet boundary conditions, they consider 
a modified Selberg Zeta function, Z (s), which has a Dirichlet series representation which 
seems to be conditionally convergent up to 5Ri = 5. For convex, co-compact Schottky 
groups Guillope, Lin and Zworski [9 1 presented numerical results for Z (5) in a large range 
of Si. They use a method based on transfer operators and due to the co-compactness they 
are able to evaluate the related Fredholm determinants in a more or less straight-forward 
manner in terms of fixed points of the corresponding maps (cf. e.g. also Jenkinson and 
Pollicott |16|). 

In this paper, we also consider an approach to the Selberg zeta function using a trans- 
fer operators. This method is applied to the family of Fuchsian groups known as Hecke 
triangle groups, generalizing the modular group. These groups have finite area but are 
not co-compact, so the evaluation of the corresponding Fredholm determinants is more 
involved. 

It can not be stressed too much that at least one of the steps in our proposed method is 
not entirely rigorous but rather supported by heuristic arguments and the entire method is 
supported by numerical evidence. 



2. Hyperbolic geometry and Hecke surfaces 

Let Jif = {z S C 1 3z > 0} be the hyperbolic upper half -plane together with the metric 
given by ds — the group of isometrics of Jf is PSL2{M.) = SL2{M.)/ {^h} where 
5L2(M) is the group of 2 x 2 real matrices with determinant 1 and I2 is the 2x2 identity 
matrix. Elements of PSL2{M.) acts on Jf via Mobius transformations. If g = (" ^) G 
PSL2{R) then z |±| and we say that we say that g is elliptic, hyperbolic or parabolic 
depending on whether | Tr g| = jo + < 2, > 2 or = 2. The same notation applies for 
fixed points of g. A parabolic fixed point is a degenerate fixed point, belongs to dJ^f and 
is usually called a cusp. Elliptic points z appear in pairs, one belongs to Jf and the other 
one is in the lower half-plane J^/f and its stabilizer subgroup in F is cyclic of finite order 
m. Hyperbolic fixed points appear also in pairs with x,x* G dJ^, where x* is said to be the 
conjugate point of x. A geodesies 7 on J^f is either a half-circle orthogonal to M or a line 
parallel to the imaginary axis and the endpoints of 7 are denoted by 7± G dJ^f. 

Let 7Z : M' — > ^ — V\M' be the natural projection map, i.e. n (z) = Fz then 7* = (7) 
is a closed geodesic on if and only if each 76 TT^' (7*) has endpoints which are conju- 
gate hyperbolic fixed points. This gives a one-to-one correspondence between hyperbolic 
conjugacy classes inF, i.e. flie set {[P] |f € F, | Trf | > 2} where \P\ = {APA-^\A e F}. 
It is known that any hyperbolic element P can be written as a power of a primitive hyper- 
bolic element, Pq, i.e. P = P™ for some m> 1. We denote this integer by m (P). In terms 
of closed geodesies on ^ this means that every closed geodesic has a minimal length 
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obtained by traversing it once only. We can now define the Selberg zeta function for F as 

(1) zr{s)^ n n(i-^(^o)-'- 

where is the set of primitive hyperbolic conjugacy classes in F, Pq is a representative in 

this class with Tr Pq > 2 and the norm of P, .yV (P) is the solution of | Tr P| = ,/K 2 + ^ - 2 
with 1 < -yV < oo. We observe that -yV (Pq ) = ^ (^o)" and since the trace is invariant 
under conjugation yV is constant over conjugacy classes. If yis the geodesic corresponding 
to P then the length of 7 is / (7) — \n.yV (P). For 5Ri > 1 the logarithm of Zp {s) can be 
written 



(2) -lnZr(i) = -^ £ ln(l-^(Po) 

= I L--^(^o)-™- — ^ 



[Po]^^ «M « 1 - ^ (PS) - ' [Pi'Sr, 'nip)i-^{pr' 

For an integer q>3 the Hecke triangle group Gq is generated by the maps S:z^ — ^ and 
T : z+ where Xg = 2cos ( - ) e [1,2). Let 7^, = — § • One can show (cf. e.g. lfT9l 



VII]) that Gq is a Fuchsian group (discrete subgroup of PSL2{M)) with the only relations 
52 = {STy = Id and which has = {z G 1 e |z| > 1 } as a closed fundamental 
domain (with sides properly pair-wise identified). It follows that Gq is co-finite, meaning 
that the Hecke triangle surface, y^q = Gq\Jf, has finite hyperbolic area. In the following 
we usually write A for Xq, Jfq and denotes the set of hyperbolic respectively primitive 
hyperbolic conjugacy classes in Gq. 

3. Symbolic Coding 

In 129)1 we showed that the geodesic flow on the unit tangent bundle of ./#^, = 
y^q X can be coded in terms of regular A-fractions (nearest A-multiple continued frac- 
tions). If {ji:};^ ^ Lx + ^ nearest A-multiple function we define Fq : Iq Iq by 
Fq (0) = and Fq{x) = -\~ {x})^ X for x ^ 0. For any number x G M we obtain the regular 
A -fraction of X, Cq {x) — [ao;ai ,a2i • • •] by first setting = ' —x — aoX, and then re- 
cursively set a„ = {Sxn]x andx„+i = Fq (x„) forn > 1. Note thatx = lim„^oo T"°ST"' ■ ■ ■ ST"" 
If X is a cusp of Gq this algorithm terminates and we get a finite A -fraction and if x is a 
hyperbolic fixed point of Gq then it has an eventually periodic A -fraction. It follows that 
Fq acts as a left shift map on £/q, the set of regular A-fractions viewed as a subset of Z^. 
If flo = we usually omit the leading "ao;", repetitions in the A-fraction are denoted by 
powers and infinite repetitions by overlines. 

In II22I it was shown that Fq is almost orbit equivalent to Gq, that is, two points x,^ G R 
are equivalent under the action of Gq if and only if, either they have regular A-fractions 
with the same tail or x has the same tail as r and y the same tail as — r (or vice versa). 
Here r G is a special hyperbolic point which can be given either in terms of its regular 
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A-fraction or explicitly. For even q one has r = \—X and (r) — [l'' ',2] with h = 
and for odd q one has r = R — X where R is the positive solution ofR^ + {2 — X)R—\ = 1 
and Zc, (r) = [l^XT^^^] with li=^. 

Let be the set of all purely periodic regular G^-inequivalent A -fractions and set 
/^r _ t^\^ {—r}, i.e. the set of purely periodic regular A-fractions with tail not equivalent 
to — r. Let £P[ denote the subset with minimal period A: > 1 and set = Ui>i 

It was also shown in |29| that for the part of the geodesic flow not disappearing into the 
cusp there exists a cross section E and a first return map =5?^ : E ^ E which has as a factor 
map in the expanding direction the map .%'.Iq^ Iq given by powers of the generating map 
of the nearest A -multiple fractions, F^. Closed geodesies on ./#^ correspond to the orbits 
of fixed points of i7 and it is easy to verify that these correspond in fact to fixed points of 
Fq, i.e. points with purely periodic regular A-fractions. It follows that there is a one-to-one 
correspondence between and -y^^. 

In practice, if Cq {x) — [oi, . . . ,«„] then = ST"^ ■ ■ ■ ST"" e Gq is hyperbolic, has at- 
tractive fixed-point X, repelling fixed-point jc* and the geodesic 7(x,x*) is closed. Fur- 
thermore y = has dual regular A-fraction (cf. f29l) c* [y) ~ [ and 
y G [-R,-r\?.gn{x). 

This connection (coding) between primitive hyperboUc conjugacy classes and periodic 
orbits of Fq is precisely what allows us to relate the Fredholm determinant of the transfer 
operator for Fq to the Selberg zeta function. 



4. The transfer operator 

In this section we will construct the so-called transfer operator for the map Fq defined 
in the previous section. 

4. 1 . Markov partitions. There is a particular Markov partition of Iq with respect to Fq 
which is important here, namely the one determined by the orbit ("'" ^) ' ■■' ^ ^ ■ ■ ■ ' "^j 
of the endpoints ±^ under Fq. Let { J^j}^.^^ be the decomposition of — ^, ^ deter- 
mined by this orbit with _/k = { 1 , 2, . . . , JC, — ff, . . . , — 2, — 1 } , J^/ = [(^j- 1 , (pj) = —-^-j 
where the order of G = |^v/^o| = {0./} J=o as - | = < < • • • < ^ic = 

and = -^j. If q is even K-=^=/!and-|== [l''] . If ^ is odd fc = ^ = 2/i + 1 
and — ^ = [l'',2, l''] . It is easy to verify that the closure of the intervals, j} is indeed a 
Markov partition of Iq for Fq. Let (p„ (y) = ST"y = then the most important property 

of the partition { } is the fact that if y e then Fq~^ (y) — {% (y) \ n e ^j} where ^ 
is a fixed set of integers depending only on j. It is now easy to show that if / > 1, / G 
and y G then 



(^r')'w= U U 



■°(Pn2° <P«i (y) 



(«!,...,«/) e ^i'j 



where we define .yVj' := (ni,n2i ■ ■ ■ ,«/) <= ^' 



ST"'ST"i-^ ■ ■ ■ ST"i Ji dJj). Let Z> 



e ^ I > m} and for A C Z let -A = e Z I - 7 e A}. It is shown in i23J that jYij e 
{{1} , {2} ,Z>2,Z>3} for /,y e ^k:, y > and that J/i^j = -jY-ij. It is also shown that 
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the non-empty elements of (for j > 0) are given by the following expressions; 

^l.lh ^ {2},^,2A+l =^>3,^-l,2h = ■yy-l.lh = {1}, 

= ^-i,2h = m^'A.lh+X = ■y^-i.lh+l = ^>2, 3 < / < 2/! + 1, 

if q is odd and 

'^^iji = Z>2, '-y^-iji = ^>i, 
^i-i - {1},^^/, - Z>2, ^-,v, = Z>i, 2 < / < 
if ^ is even. For example 

(-^.) - (S :£) ' foi-'? = 3. (.y^;) = (g :g) for^ 4. 

4.1.1. Transfer Operator corresponding to F^. For any interval / C M let C (/) denote the 
space of continuous real- valued functions on /. If / G C(lq) the transfer, or generalized 
Perron-Frobenius operator corresponding to Fq, acts for real j3 > 5 on / by 



where is the characteristic function of J^j. It is important to note here, that ££^f [x) 
is in general not continuous, but only piece-wise continuous. For this reason consider 
the action of on vector-valued functions in = 0;g^^C(^). For f G we set 

/ {^) = fi {^) if G ^i- We can now write : as 

(^pf),{x) = j: j: \K{x/fj{(Pn{x)),iej, 

^ ' ne.A<j 

respectively, for any I > 1 
where 

^^..fj (x) = K^^"' • • • ST"')'xffj (5r'i • • • ST"'x) . 

(ni,..,«,)G../t;j 

It is convenient to use a composition operator Ttp related to the principal series repre- 
sentation of PSL2iRy. Define 7tp{A)f{x) \A'{x)f f{Ax) = \cx + d\'^l^ f {!^) for 
A e P5L2(IIi)- Note that tt^ {AB) = Ttp (B) Ttp (A). With this notation one gets 

(«l,...,n,)G./^j 

To obtain better spectral properties for the operator we have to restrict its domain of 
definition even more. For any open disk D in C we let M {D) be the Banach space of func- 
tions holomorphic in D and continuous on the closure D together with the supremum norm. 
Let {Di}-^j!^ be a set of open disks with diameter which contains an e-neighborhood of 

J^i, constructed in such a way that for (ni ,...,«/) G one has o (p„^ ° - ■■°(Pn, (A) C 
Dj. That such a choice is possible is shown in ll23l . Let = ^ (D,) and define the Banach 
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space !M = ®iej with norm given by 



maxy ||/; II Then we want to consider 



as acting : ^ ^ For this purpose we also need an analytic extension of . 
If A = ^) e P5L2(K) and A{D)<ZD for some disk D then TZp (A) ■.S§{D)^S§ (D) is 

defined for any j3 e Cby iA)f{z) = (^{cz + dy^Y f (^^). Usually we simply write 

the first factor as {cz + dy^^ , but remember that there is a choice of sign involved, i.e. 

{-cz~dy^^ = {cz + dy^^. For/> 1 and/G ^: 

(3) i^f) .{z)= ^ ijfj iz) , i e A with 

We now have a representation of the operator asa(K"+l)x(K'+l) matrix of operators 

Next we need some facts from Grothendieck's theory of Fredholm determinants and 
nuclear operators on Banach spaces |8| (Ruelle 1261 provides more detailed references). 
The following Lemmas follow from this theory. 

Lemma 1. Let D be any open disk in C and let ^ (D) be as above. IfW : ^ {D) — > ^ (D) 
is a simple composition operator^'f [z) = \if {z) f {z)) with \j/,(p continuous in D and 
(p (D) C D. Then (p has an attractive fixed-point z« S A f is nuclear of order zero and has 
trace Tr pj(o)'i' = T^^jh)- 

The formula for the trace, sometimes referred to as a special case of the Atiyah-Bott 
trace formula is easy to verify directly since the eigenvalues of ^ are all of the form /i„ = 

V/(z*)(9'(z*))",«>0and |^'(zh.)| < 1. 

Lemma 2. If ^ is a nuclear operator of order zero on a Banach space we can express the 
Fredholm determinant det (1 — ^) in two different ways: 

- logdet (1 - ^) = £ i rr if ' = - log n (1 - Xj) 

1=1 ' ,;=1 

where {Aj} J_j are the eigenvalues of^ (countedwith multiplicity). Furthermore, if^ = 
if (s) is a meromorphic function of s then det(l — i^(s)) is also meromorphic in s. 

Proof Cf. e.g. 1 8 , prop. 1, pp. 346-347]. □ 

Lemma 3. Let A G SL2{R) be hyperbolic with attractive and repelling fixed points xj^ and 
x^ respectively. IfD is a disk with diameter on R containing only the attractive fixed point 
of A then A{D) C D. 

Proof. This Lemma is easy to verify by conjugating with the map in 5'L2(IIi) which takes 
Xj^ to 0, x_ to i°° and A to z i-^ fiz with < / < 1 . □ 

If A = ( " J ) is hyperbolic with attractive fixed point x+ it is easy to verify that ,yK (A) = 
jA where ja {x) = cx + d. Since TZp (A)/(jc) = Ja {xy^^ f (Ax) it is easy to see that 
ifx+^D and x^ ^ D then by Lemma[T]7r|g (A) is nuclear of order zero and 

.yK{Ayl^ 

Tr^(.)^/3(A) = ^-;^. 
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Let « = . ,«,) e and set ^ ST"^ST"^ ■ ■ ■ ST"'. ThenA,iJ^; C J^j so the attrac- 



tive fixed point of A3, Xj^= [«i , . . . , «/] G and by Lemma[3]it is clear that A^ (D^) C Dj. 
This demonstrates that all composition operators showing up in the operators Jf^ jj ap- 
pearing in the trace of are nuclear of order zero for 5Rj3 > 5. The arguments in ^24\ or 
lITSi can be generalized to show that is also of trace class and nuclear of order zero for 

9^/3 > J. It is clear that Tr ^^p = Y^ie/x .^/-^^ /, ^'^^ ^ — ^ ^^^^ ^'^^ '^^'^'^ 

to see that any off-diagonal term, .j : S^i is a bounded operator. One can now 

use similar arguments as those in Q to show the following lemma. 

Lemma 4. If SRj3 > | then is nuclear of order zero and hence of trace class. 

From the identification of hyperbohc conjugacy classes with purely periodic A-fractions 
we may now calculate the trace of I£p 

Tr^^^ = £Tr^,^.|.^=£ £ Tr ,^,7r^ (A„-) 
y .yK{ST"^ ■■■ST"'y' 



«l,n2, ■■■,«/ £.5*^; 



l-.yK(5r«i •••5r"0" 



By the standard Grothendieck theory the Fredholm determinant of 1 — S^p is well-defined 
and can be calculated by 



/=! ' Z=I 



We now need to study the relation between A-fractions and hyperbolic conjugacy classes 
in more detail. Let x = [ni, . . . ,«/] correspond to the hyperbolic A=Aj, = ST"^ ■ ■ ■ ST"' = 
Pq where Pq is a primitive hyperbolic and m = m (A) then x = n\,. . .,nig where Iq is 
the minimal period and Iqiu = I. Furthermore, all shifts, [«/,n,+i, . . . ,n,_i] for 1 < / < 
«/p belong to the same conjugacy class [A] and the norm is also constant over conjugacy 
classes. Let ^^{x) — ^ {Dj) where x S Dj, set A^r = r and ,J(f^ — (Ar). Then 

- logdet (1 - .^p) +logdet (l - Jfp) 
= I y L Tr (Ah) - £ j Tr ,^(,) np {a[ 

= lj E Tr,^^.^^,;r^(A„-) 

nGZ',.t=[«i,n2,...,H;lG.5J"" 



'0 = 1 )1o£Z'o^x=[ »i ni^ ]e.9"' ^ ' 

= E^ E JkT^'^^x«o)^/^(^") 

/o=l'0[Po]e.ir,''^^ 
^ 1 .-^(P)-^- 



m(P) 
no 



g 
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and by compaiing with ^ we see that for 9^/3 > 1 we have \nZg ()3) = logdet (l — Jf|3) 
logdet (l — Jtp^ and thus 

det(l-^«) 

(4) ZJp)^ q-. 

detfl -jr«) 



But since the right hand side is in fact meromorphic for j3 e C this equation provides an 
analytic continuation of Z (j3 ) for j3 e C. Note that the operator Jtf is a simple composition 
operator and we can evaluate det (l — Jtj^^ explicitly. It is easy to show that all eigenvalues 
of J^^ ai-e of the form fi„ = {2 + A/;)"^'"+^' , n > and hence 

(5) det{l-Jfp) = Yl{l-^„). 

ii>0 

To evaluate the factor det ( 1 — J^'p ) we use the same identity det ( 1 — ) ~I\n>i{^ ^ ^n) 
where {^n}„>i are the eigenvalues of J^p counted with multiplicity. In the next chapter we 
will discuss how to calculate eigenvalues of ^p . 

Remark 1. The method to relate Zp {s) to Fredholm determinants of nuclear operators can 
be extended to any finite dimensional representation % of T. The identity (01 will hold 

with Zr, ^p and Jfp replaced by Z^ {s) = W^p^} Ilirdet - j [P) .vV {Py''^ , respec- 
tively and J(fp. Here if^ and Jf^ ai-e obtained by replacing np (A) with (A) ^ 
X (A) 7tp (A) in all formulas. The only problem is to obtain explicit expressions for the 

truncated operator si^p^^ which will be introduced in the next section. The algorithm has 
been implemented and tested for representations induced by the trivial representation of 
the Hecke congruence subgroups Fq (p) with prime p. This allowed us to compute e.g. 
2ro(p) (i)for/:' = 2,5. 



5. Analytic Continuation of and Computation 

It turns out that the same analysis which enables us to deduce an analytic continuation 
of ^p to j3 eC is also vital to compute the eigenvalues of ^p. 

We follow the same procedure as in e.g. ("ZS' Tl to demonstrate that Jiffp admits a mero- 
morphic extension to the whole complex plane. First of all we have to change domains 
once more. To make some of the calculations easier it is desirable to work with functions 
which have power series expansions around zero. For this purpose we choose open disks 
Dj D Dj such that G 5, and ST"Dj C Dj for n e ^ij. That this construction is possible is 
shown in ||23l . If = 3§ (5,) then ff G !Mi has a power series expansion centered at zero. 

If / e it can be shown that either .J/ij = {riij e Z\ {0}} or = for 1 < ; < 
fc — 1, -yViK = {« G Z I n > n,K-} for some > 1 and ^i-j = —jV-ij. Let 1 < / < fc and 
consider .^p^ij ■ ^i- Let / G and > 1. Taylor's theorem with remainder gives 
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/ {z) = iLo akZ^ + (z) with Rn {z) = 0[ \z 

^P,iKf{z) 



N+\ 



Then 



ne^K n>niK 
2p 



I 



n>nu 



1 



z + nX 
1 

z + nX 



f 



-1 



z+nX J 



= ^^lf{z)+^'p;'j{z 



2/3 



£ayi 



-1 



z + nX 



+Rn 



1 



where Sf^'^^f (z) = [/ (z) - if^i aA:Z^] is analytic for 9^j3 > and in fact nuclear 



of order 0. We also have 

AN) 



operator on the other hand can be written as 



N 



j_|29?/3+iV+l 
nX I 



as AT ^ oo. The 



^z + nA , 



where ^ (i,z) is the Hurwitz zeta function. It is known that for any z S C the function 
C, {s,z) is meromorphic with only one simple pole at 5 = 1 with residue 1. Hence sz^^i^ is 
of finite rank and meromorphic in j3 with at most simple poles at the points /J^ = =^ 
0<k<N. For the operator corresponding to c/ff,- = {tiij} one has 



2 ' 



7tp{ST'"J)f{z) = {z + mjX)-'^f 

N / 1 \ k 



-1 



= {z+HijXy 



-2J3 



k=0 



-1 



z + riiiX 



z + riijX 
+Rn 



1 



z + riifX 



p,'j 



fiz) + ^^p{z) 



where s^^'^^fiz) =lto«* (-1)' iz + nur'^'P ^d^^^.fiz) = {z + mjX)-'pRN 
o(\z + nij 



It is clear that in this case , 



. ' is entire of finite rank and that 



is entire and nuclear of order 0. 
Since N>1 was arbitrary, we conclude that all components ^pjh have meromorphic 
continuations to the entire complex plane with at most simple poles at the points Pk = =^ , 

k = 1, . . .. The same clearly holds true for the operator ^p. Note, that in the determinant 
det (1 — poles may well cancel against zeros due to the presence of eigenvalues equal 
to one. 



5.1. Computation of ^Z^. Let = s^^ 



AN) 



components -^p^jl respectively -^^'^^ given above. To obtain a numerical approximation 
of det (l — Jfp) it is necessary to approximate the spectrum of ^p. For this purpose we 



AN) 



-^^p where j^^^ and have the 
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construct another finite rank approximation of -f^p^^ in terms of a matrix which is more 
suitable for computations. 

Let be the space of polynomials of degree less than or equal to A^. Then is a 
subspace of all ^,'s and we let 11^ denote the projection given by truncation 

of the power series, i.e. n^v (£"=o^fc^^) ~ £^o^fe^^- ^'^ denote 

the projection from SS = ©,e^j. to the space ©ie^^. obtained by truncating each 
component. 

We saw, that ^^^^ maps ©,g^^ —■ into a space spanned by Hurwitz zeta 
functions. By applying Iliv to the resulting expression we obtain an operator ^/^^'^^ : 

£2k(n+i) _^ £2k{n+i) yyijich can be represented by a jcat x Kjy complex matrix where Kjy = 
2k' (A^ + 1 ). This construction will now be explained in detail. 
Let > 1 be a fixed integer, then with the notation as above 

N 

n=0 k=0 

where Ui^^nk = oiiK,nk iP) = „|^„+^+2/j (^ + 2j3)„ C (2j3 +k + n,niK). For i e and 1 < 
7 < ff — 1 we get 



k=0 



k=o n=o n\ {nijA) '^^ ^ 



oo N 

n=0 k=0 



where a^-^j: = ay,„fc (J3) = ^[^^z^j^t+n (2j3 + A:)„ n,/^ * ". For n < -1 we use a slightly 
modified definition of TT/j , namely {ST-") f (z) = + 



It is then easy to see that 0£,;,„ft = (—1 )"''"* cti-j^nk for 1 < 7 < k"- 

By truncating the sum over n at A'^ in the formula for s^^ij we get operators -^pl^ '■ 
- and ^j'^'^) = {^^"^if^) . .^^ . Then ^^''^ = ^yj ° Oiv and ^^(^''"^ = o 

o Un and by the identification ^ €^'^(^+1) it is clear that ^(^-'V) : ^ -^w can 
be represented by the Kjy x jqv complex matrix 

= {"ij,nk) ijejr^,0<n,k<N ■ 

In the next section we wiU discuss the relation between the eigenvalues of ^ and those of 



COMPUTATION OF SELBERG ZETA FUNCTIONS ON HECKE TRIANGLE GROUPS 



II 



5.2. Approximation of the spectrum of LetE|3 denote the spectrum of the operator 

and the spectrum of si/^^''^^ . Since is nuclear, setting II'^ = Id — 11^ then 
is bounded and it is easy to verify, that the conditions of Theorem 2 and Proposition 3 in 
Baladi and Holschneider [3 1 are satisfied by the approximations si^p^''^^ ,N>1. Hence the 
following Lemma can be deduced: 

Lemma 5. Let f ^ be an eigenfunction of Jifp corresponding to the eigenvalue Xp G "Lp 
with algebraic multiplicity d. Then there exists Nq >0 such that for all N > No there exist 
eigenvalues X^j G T/p with corresponding eigenfunctions f^j , I < j < I such that the sum 
of the algebraic multiplicities of Xnj equals d and 

max ( I A - Xnj \ , f-fy.] ) <c{N), 
i<j<i\' ' / 

where c {N) as N °°. 

Definition 1. If Xnj G E^ is one of the eigenvalues in Lemma|5]approximating a Xp ET.p 
then Xfi.j is said to be regular, otherwise it is said to be spurious. 

A problem in computing the spectra of ^p using s^p^'^^ is that we do not know a priori 

which eigenvalues of s/p^'^^ are regular and which are spurious. A trivial consequence 
of Lemma|5]is the following Lemma which gives a necessary condition for a sequence of 
eigenvalues X^.j. G T!p to be regular. 

Lemma 6. Let {-^a^, j;},>i be a sequence of eigenvalues of s^p^'''^'^ such that Xn^j- — > A G 
T.p as j — > CO. Then for any e > there exists No>0 such that \XNijj — XmiJi \ < £ for all 
Ni,Mi>NQ. 

Remark 2. Let T.p = {Xp „}^^^ (where eigenvalues are counted with multiplicity). By 
Bandtlow-Jenkinson mis] there exist positive constants A, c such that \Xp„\ <Ae^™. Nu- 
merically we fond that a similar bound seems to hold for the operators s/^'^^ . It follows 

that is a limit point of E|3 and there exist many sequences of spurious eigenvalues {XNj.jj } 
converging to 0. 

We wiU now present an algorithm which uses Lemma|6]to compute an approximation to 
[s), but to put Lemma|6]into praxis we first need to make the following heuristic claims. 

Claim 1. There is no sequence Xn^j-, G > 1 such that Xn-^j^ — > X unless A G Ejg or 
A =0. 

Claim 2. Suppose that = {■^a',(}i<,<ic^. |Aa?,„, | > |AAr,„2 1 > • • • > IAat.h^ | > are regular 
eigenvalues with an estimated error < e and |AAf,„j^| < e. Then there does not exist an 
eigenvalue Xp G E^ in the region {z G C| |z| > |AAr_„^| , max,=i_..._/f|z — AAr,„,| > e}. I.e. 
the sequence {XN,ni}\<i<K approximates all eigenvalues of J^p with absolute value greater 
than or equal to \XN,nf, \ . 

Algoritlim. Let 5, e > and consider N and M for some M>N+\. 

Step 1: Compute the two spectra E^ — {■^Af,i}i<,<Kjv '^"'^ ^ {^^'j^ KjKkn (both 
ordered with non-increasing magnitude and repeated according to multiplicity) and the 
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Step 2: Letk — and consider in sequence each i — I, . . . ,Kfj. If there exists a j such that 
5i j < 5 we assume that X^jf. and ^m./V isfe approximating some A G E^g and accordingly 
increase k by 1, set ij^ = /, ji^ = j, Si^ = 5;^^^, and Xp j. — Xmj,^- 

Step 3: Let K denote the last value of k. Then {Apj-j^^j is an ordered sequence of 
eigenvalues believed to approximate eigenvalues ofJifp and we define 

^V,M(i3) = n(l-^i3,A)- 

.7=1 

If > Ewe increase N andM and start from Step 1. As will be explained in section W2\ 
below it might also be necessary to increase the working precision simultaneously with N 
in this step. If \Xfj [^\ < £ we assume that d^.M approximates det (l — and return 

ZgiP)^ Zg,N.M iP ) - dNM (i3 ) det ( 1 - ) " ' 

as a tentative value ofZq (s) with an assumed error depending only on 5,£ and the working 
precision. The factor det(l — J^p) can be computed using relation Q to any desired 
accuracy. 



6. Discussion of the Results 



The numerical method. Algorithm 15.21 which is proposed as a means to evaluate the 
Selberg zeta function relies on the heuristic Claims [T] and |2] above. It is thus clear that 
no amount of internal "consistency tests", e.g. stability under change of order of approxi- 
mation and variation of the parameters e and 5, can certify that the result returned by the 
algorithm is correct. If Claim[T]is wrong we would obtain extra eigenvalues not associated 
to ^p and on the other hand, if Claim|2]is incorrect we might actually miss eigenvalues of 
comparatively large magnitude. In both cases we would only be able to approximate Zp (s) 
times some unknown factor. 

The need of an independent test to verify the accuracy of our numerical results is thus 
obvious. We propose to use a test relying on the functional equation of (s). The setup 
will be discussed in Subsectionl6.1l 



Remark 3. If we were only concerned about zeros of Zp on the real axis, i.e. eigenvalues 
equal to 1 of ^p for real j3 much more is known about approximation of eigenvalues and 
eigenfunctions, cf. e.g. lfT7l(20l . 

6.1. The functional equation for Z^ (s). Let Zg (s) be the Selberg Zeta function for Gg. 
We know HIl p. 499] that 

(7) ^^^^,{s)cW,{s), 

where (pq [s) is the scattering matrix (here a 1 x 1 -matrix), c = = ±1 and 



>P (i) = — — ^exp n ttw{nt)dt+ 

r(i + i) \ q Jo 



ci-l , ... 1 



^^^t'l^sin^io [Y+T^^Y+e^ J ' 
(1 -2i)ln2). 
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The function (s) can be computed to any desired degree of accuracy using standard 
methods of numerical (e.g. Gauss) quadrature. Evaluation of (Pg (s) on the other hand is 
more tricky. For q ~3 there is an explicit formula 1121 p. 508] 



r(,v-i)C(2.-i) 



(8) (p^is)^^- ^,.f,i^. 

T{s)C,{2s) 

For q>A the only explicit formula is in terms of a Dirichlet series with abscissa of absolute 
convergence equal to 1, cf. e.g. ifTZl p. 569] and ll30l . Note that for <7 = 4,6 it might still 
be possible to work out explicit formulas for <Pq{s) using the relations between G4, Gg 
and Fo (2) , Fq (3) respectively. We do not pursue this approach and for all ^ > 4 we use 
values of (pq (s) obtained by an algorithm of Helen Avelin [jT]. The main idea of Avelin's 
algorithm is that <Pq occurs in the zeroth Fourier coefficient of the Eisenstein series E (s,z) 
for the group Gq and one can use a method based on the G^— invariance of E{s,z) to 
compute its Fourier coefficients and thus also (j?^ (i). This method was first introduced to 
compute cuspidal Maass waveforms on Hecke triangle groups by Hejhal fTSl [T4 1 . Later it 
was generalized to the setting of general subgroups of PSl^CZ) (2^, Ch. 1] and finally it 
was generalized to compute Eisenstein series on Fuchsian groups with one cusp by Avelin 

m. 

Another application of the functional equation is that we may define a real- valued func- 
tion 

(9) ^qit)^Zq(^+it)e-'®^'^ 



where © (f ) = ^arg ( | + /f) ^ -I- /f )) and the branch of the argument is chosen so 
that 3fq (/) becomes continuous. Note that a single choice of a branch cut is in general not 
possible because (p (^j + it^ ^ (^j + it^ winds around zero as f G varies. The advantage 
of considering 3fq (f) is in our case purely aesthetic, in that we may plot graphs of {s). 

It is known lfT2l p. 498] that Zq (s) is zero for s = si^ = j + ir/, where \+ is an 
eigenvalue of A and at ^ = 1 — 7 where <Pq (7) — 0. In Figures [T][3] we plot 2fq (f ) together 
with blue vertical lines at t = ri^ and green vertical lines at f = 87. The zeros of Zq (s) on 
and off the half-line are clearly visible as zeros and "dips" of 3fq{t) at the corresponding 
points. The eigenvalues of A were computed by the method of Hejhal indicated above, see 
e.g. 11131 [1411281 and zeros of (pq (s) were located using Avelin's algorithm. 

Verification of these zeros as well as the zeros on the real axis of Zq (s) ( lfT2l p. 498]) 
does of course also lend credibility to our proposed algorithm but since this verification 
does not tell us anything about the accuracy for general s we prefer to concentrate on the 
error estimate using (pq (s). 

Remark 4. The actual value of (Pq{j) G {±1} can be computed experimentally in two 
different ways. The straight-forward way is to use Avelin's method but it is also known 
(cf. (IT, p. 498]) that Zq (5) has a simple pole at 5 = j if and only if (pq ( = — 1. Experi- 
ments performed using both methods indicate that (Pq{j) = — 1 for all q' > 3. 

In certain cases one can use the transfer operator to show that Zq (s) has a singularity 
at s = I by showing that J^fp has an eigenvalue jlp ~ i| in a neighborhood of p — j, 

but it is not possible to exclude that this singularity in det(l — ^/s) is canceled by the 

appearance of an eigenvalue = 1 for Jffi . 

1 

6.2. Discussion of data and error analysis. The procedure for testing and producing 
error estimates of the proposed algorithm to compute Zq (s) is now clear. Given tentative 
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values of Zq {s) and (1 — s), denoted by Zq (s) and Zq (i) we compute the quantity 

%{s)^Zq(\-s)Z,{s)-'^{s)-' 

and compare this with the value of c (pq (s) obtained as described above (in all cases con- 
sidered here we have c = 1). The difference \(pq (s) — ^q (s) | or relative difference in the 
neighborhood of a zero of (p (s) gives an estimate of the accuracy of the values Zq (s) and 

Zq{l-S). 

To confirm a value Zq (s) we thus need also to compute Zq{l ~ s) , but on the critical 
line with s — j + it we have Zq{ \ ~ s) — Zq {s) so we need only compute Zq [s). 

Using this "^-test" we may verify the correctness of Claims [T] and |2] As it turns out, 
these two claims seems to be correct in theory. In practice, however, they and the entire 
algorithm may fail unless the working precision is increased as necessary. This phenome- 
non is clearly visible in Table [T] where we investigate the case q — i and s ~ ^ +5i using 
different degrees of approximation (here M =N + 3 always) and working precision WP. 
In this table we list the estimated error, | (p^ (s) — q)^ (s) | , the time it took to compute Z3 (s) 
in seconds, the number of eigenvalues of Jfj3 which were used in the computation, the size 
of the smallest of those eigenvalues and the maximum of differences between eigenvalues 
of and With working precision of 50 digits we see that the error decreases 

as N = 25,50 and 75. To further increase up to 100 does not improve the accuracy and 
increasing up to 200 actually results in a worse approximation than at = 25. The 
reason for this phenomenon is that Claim [T] is violated due to an increasing number of 
spurious eigenvalues and in particular there appear spurious eigenvalues which do not vary 
fast with A^. This problem can be overcome by increasing the working precision, which 
is demonstrated in the remainder of the table, where the precision has been increased to 
WP = 100, 150 and 200 digits respectively. To know a priori when the precision has to be 
increased one must study more closely the relative differences 5k. For example, in the case 
WP = 50 and A^ = 200, the relative differences for the spurious eigenvalues of a certain 
magnitude are much larger than the relative differences of regular eigenvalues. If one sees 
such a break from the otherwise almost monotonously increasing 5it it is a clear sign to in- 
crease the working precision. What is not visible in this table, is that the need for increase 
in precision is actually dependent on the matrix size Kjv = 2 K" (A^ + 1 ) and not only on A^. 

Table|2]contains values of 93 (5 + «/) , 1 < n < 10, computed using A^ = 100, M — 103, 
5 = 10^^ and 100 digits working precision. The third column contains the true error, i.e. 
(pi compared to the explicit formula ([8]) for ^3. One can see that in this case the true error 
agrees well with the error estimate in the fourth column given by the absolute value of 
the smallest eigenvalue used in computing d {s). The fifth column contains the difference 
between (pj (s) computed by Avelin's method (using double precision) and by the explicit 
formula. 

In Table 16.31 we list values (p^ (s) of (p^ +«/) , for 1 < n < 10 given by Avelin's 
algorithm and of ^4 {s) by our algorithm using A^ = 100, M = 103, 5 = 10"^ and 100 
digits precision. The fourth column contains the difference between these values and the 
fifth column contains the size of the last eigenvalue Xk used in the evaluation of d^M 



Comparing the values of 



Xk 



in Tables |2] and 16.31 we observe that the eigenvalues of 
for q = 3 seem to decay more rapidly than for q = 4. We would also expect the errors 
in the tabulated approximations of ^4 (s) to be greater than those of ^3 {s) even though we 
use the same level of precision and approximation. However, there is no reason to believe 
that the error in (p^ is any worse than in (p^. Moreover it is very unlikely that the values 
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(p^ (s) should agree with our values ^4 (s) to a much larger degree than the true accuracy, 
cf. e.g. s = i+ 9/ where \1k\ =2-10-^'' but | (P4 W - ^4 (^) | =2 - 10"'^ We conclude 



that the values of 



do not give an accurate estimate of the true magnitude of the error 
for q = 4- but that they still provide us with an upper bound. 

The final conclusion we can draw from Tables [11 16.31 is that the value of \Xk\ alone 
is not enough to estimate the error in Zg {s) unless the working precision is high enough. 
To completely eliminate the external test by using (pg to confirm values produced by our 
algorithm one needs a better understanding of when it is necessary to increase the working 
precision. 

It is clear, that high precision eigenvalue computations are very time consuming. To 
evaluate Zg (s) to a fixed precision it is necessary to increase the approximation level 
as 3 s grows and this forces a simultaneous increase in the working precision. Altogether 
this makes it very time consuming to compute values of Zg + 't) for large f's and to 
reach even values of f « 1000 for q — 3 seem to be out of reach with current methods and 
hardware. Remember that the size of the matrix Kj^ grows with q, so similar problems arise 
when computing Zg (s) for large <7's. To end this discussion I would like to give a feeling 
of the necessary CPU-times. 

To compute Zg (5 + it) with an estimated error of 10^^ for q — 3 takes 12 seconds at 
f = 15 and 30 minutes at f = 100. Decreasing the error to 10^^ at f = 15 only increases 
the time to 26 seconds. For ^ = 8, to compute Zg (5 + 15/) with an estimated error of 
10~^ takes 5 hours and 57 minutes. If these figures seem outrageous, remember that for 
<7 = 8 the size of | = 6 and in this case = 190 so K"Af = 1 146 compared to Kn = 352 
for ^ = 3 at f = 100. We see that ^ w 3.25 and = « 1 1 .9 w 3.45^ so the 

CPU-time increases roughly like the square of the size of the approximating matrix, which 
is to be expected from the eigenvalue computations. 

6.3. Implementation. The algorithm outlined on p. [TT| above has been implemented in 
Fortran 90 using the ARPREC |2| library for arbitrary precision computations. It was 
also necessary to write arbitrary precision Riemann and Hurwitz zeta functions as well as 
an arbitrary precision version of the standard linear algebra system LAPACK. Fortran 90 
codes can be made available from the author upon request. 

For the interested reader who is not comfortable with Fortran there is a version of the 
algorithm implemented in MuPAD ifTSl and this version is available from the homepage 
of the author The choice of MuPAD is mostly because of its good multi-precision linear 
algebra capabilities. 

Avelin's algorithm is currently only implemented in double precision FORTRAN 77 
hence the MuPAD version only contains the complete error check using (pg (s) for q — 3. 



Table 1 . Demonstrating interplay between working precision and level of approximation in computing Z3 {s) 



WP 




Z3(|+5i) 


1^- 


<pl 


Time (s) 


K 




\Ik 




maxj^ 


5k 


50 


25 


1.1954 + 0.0811/ 


1 


• 10 


-2 


12 


3 


5 


10- 


-2 


8 


10 


-9 




50 


1.192213397499979 +0.074413721696096J 


5 


• 10 


-9 


65 


12 


2 


10- 


8 


5 


10 


-8 




75 


1 . 192213402687674941 183 + 0.07441372136992775? 


7 


•10 


-13 


202 


17 


3 


10- 


12 


2 


10 


-8 




100 


1 . 19221 3402687674941270 + 0.074413721369927750/ 


7 


• 10 


-13 


839 


19 


3 


10- 


12 


2 


10 


-8 




200 


1.0165+0.0782; 


3 


• 10 


-2 


3764 


19 


4 


10- 


-9 


5 


10 


-8 


100 


100 


1.192213402686855883038325 + 0.0744137213702737315168/ 


4 


• 10 


-19 


599 


26 


4 


10- 


18 


2 


10 


-10 




200 


1.192213402686855883047193363+ 
0.0744 1 372 1 37027373 17790 1 83373/ 


3 


•10 


-25 


4324 


34 


1 


10- 


-24 


2 


10 


-8 




250 


1.1837+0.0575/ 


3 


• 10 


-2 


10164 


27 


4 


10- 


18 


3 


10 


-8 


150 


200 


1 . 1 922 1 3402686855883047 1 9355 1 1 30 1 1 7623672+ 
0.07441372137027373177901851156938182265/ 


1 


•10 


-36 


6637 


50 


5 


10- 


-36 


1 


10 


-9 




250 


1 . 1 922 1 3402686855883047 19355 1 1 30 1 1 7623672+ 
0.07441372137027373177901851156938182265/ 


1 


• 10 


-36 


12847 


50 


5 


10- 


36 


2 


10 


-9 


200 


250 


1.192213402686855883047193551130117621955253934465021290+ 
0.0744 1 372 1 37027373 1 7790 1 85 1 1 56938 1 823 16447232 1 1 53746259/ 


6 


•10 


-51 


18788 


68 


7 


10- 


49 


8 


10 


-8 



Table 2. Comparing values of (p3{s) together with different error estimates {N = 100, WP = 100, 5 = 10 ^). 



n 


93(*) =Z3(1 -s)/Zi{s)/'^3{s), 5 = i+m 


\9i- 


93I 


K 




1^ 




max 


5k 


l<P3 - 


9i\ 


1 


0.523127151694381217718-0.852254646898521675788J 


1- 


10" 


-20 


27 


6- 


10- 


21 


8 


10 


-9 


8 


•10 


-17 


2 


0.777709870863430801402-0.628623382289893657301/ 


3- 


10- 


-21 


24 


1 • 


10- 


20 


3 


10 


-8 


1 


•10 


-16 


3 


0.810307536439650550895-0.586004860380103327530; 


1 • 


10- 


-20 


24 


9- 


10- 


21 


2 


10 


-8 


4 


•10 


-16 


4 


0.784 1 1 6026298 143660937 - 0.6206 14258056007668073; 


4- 


10 


-20 


27 


1- 


10- 


19 


9 


10 


-8 


6 


•10 


-16 


5 


0.620614258056007668073-0.709907649199078141317; 


4- 


10" 


-19 


26 


4- 


10" 


-18 


2 


10 


-10 


7 


•10 


-16 


6 


0.473769476721985155012-0.880648898782356035905; 


4- 


10- 


-20 


26 


2- 


10- 


18 


2 


10 


-9 


1 


•10 


-15 


7 


-0.982666838048427466635+0.185380380299279850669; 


2- 


10- 


-19 


26 


1 • 


10- 


18 


9 


10 


-9 


1 


•10 


-14 


8 


0.947280945444430850195 - 0.320404136049934947281; 


3- 


10- 


-19 


26 


2- 


10- 


18 


1 


10 


-8 


2 


•10 


-15 


9 


0.678702274737248216706-0.734413522660418366220; 


1 • 


10- 


-18 


26 


3- 


10- 


18 


1 


10 


-8 


1 


•10 


-15 


10 


-0.063355766687361081600-0.997991005384044865561; 


2- 


10- 


-18 


26 


5- 


10- 


18 


2 


10 


-8 


5 


•10 


-15 




Table 3. Comparing values of 94 (s) (A' = 100, WP 


= 100, 5 


= 10- 

















n 




q>4{s)=Z4{l-s)/Z4{s)/^4{s) 




^4 


K 




Ik 


maxSfc 


1 


-0.2632601 86 1 373 177 - 0.96472486979 1 8721 ; 


-0.2632601 86 1 373 176 - 0.96472486979 18723; 


2 


10 


-16 


27 


5^ 


10- 


14 


7 


10- 


-8 


2 


-0.7021440712594831-0.7120349030736887; 


-0.7021440712594827-0.7120349030736895; 


9 


10 


-16 


25 


3^ 


10- 


13 


5 


10- 


-8 


3 


-0.99 12520623526865 + 0. 1 3 198238095 1 195 1 ; 


-0.9912520623526863+0.1319823809511939; 


1 


10 


-15 


24 


6- 


10- 


12 


5 


10- 


-9 


4 


0.2148427612152942 + 0.9766486512320531; 


0.2148427612152920 + 0.9766486512320534; 


2 


10" 


-15 


26 


3- 


10- 


13 


8 


10- 


-8 


5 


-0. 8749676464424498 - 0.484 1 8 1 389232342 1 ; 


-0.8749676464424484-0.4841813892323440; 


2 


10 


-15 


24 


3- 


10- 


12 


3 


10- 


-8 


6 


-0.0732387210885128 + 0.9973144387470341; 


-0.0732387210885146+0.9973144387470377; 


2 


10 


-15 


24 


2- 


10- 


11 


1 


10- 


-8 


7 


-0.0299908075389591-0.9995501745601175; 


-0.0299908075389498 - 0.9995501745601 176; 


9 


10" 


-15 


24 


A- 


10- 


11 


3 


10- 


-8 


8 


0. 85545989 1 6720 125 + 0.5 17869070075 199 1 ; 


0.8554598916721716 + 0.5178690700748949; 


3 


10 


-13 


24 


1 • 


10- 


10 


1 


10- 


-8 


9 


0.7163471899280185-0.6977440099938026; 


0.7163471899280196-0.6977440099938012; 


2 


10 


-15 


24 


1- 


10- 


10 


7 


10 


-8 


10 


-0.7358033312663973-0.6771952877104747; 


-0.7358033312663925-0.6771952877104797; 


7 


10" 


-15 


24 


3- 


10" 


-10 


1 


10" 


-7 
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Figure 1 . ^3 (f ) 
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